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Abstract 

Resonances driven by the error field multipoles in the magnets 
often show themselves as islands in phase plots of the particle motion. 
In order to be able to measure the strength of the island resonance, 
and then to correct it, it is helpful to study how these resonances 
show themselves. In particular , one can study how the presence of 
the island resonances distort the tune dependence on the amplitude, 
and the emittance growth caused by the island resonance. Islands in 4- 
dimensional phase space, unlike islands in 2-dimensional phase space, 
are not easy to visualize. This study will show how to visualize the 
islands in 4-dimensional phase space by studying the tune dependence 
on the amplitude of the particle. These effects are to some extent 
measureable, and can lead to a way to correct the resonance. 

Introduction 

Resonances driven by the error field multipoles often show themselves as 
islands in phase plots of the particle motion. This happens when there is 
another non-linear field present, besides the error field multipoles, which 
produces a strong enough dependence of the tune on the amplitude of the 
particle motion. Non-linear fields that produce a dependence of the tune 
on the ampliude can be the quadrupole fringe field , space charge fields and 
octupole correctors. In some cases it may be advantageous to introduce a 
non-linear field that produces a dependence of the tune on the amplitude 
of the particle motion. In order to be able to measure the strength of the 
island resonance, and then to correct it, it is helpful to study how these 
resonances show themselves. In particular , one can study how the presence 
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px vs x, 3nux=1 9 resonance 
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Figure 1: p x vs. x for the 3v x = 19 resonance excited by a ran- 
dom b 2 in the SNS magnets. The initial coordinates have Xq = 24mm 
to Xq = 44mm in steps of 2 mm with p x o = and e y = .u x q = 
6.3267, Vyo = 6.2267. In the figure x0,px0, y0,py0, epxO, eptO, nuxO, nuyO 
represent x , p x0 , yo, p y o, e x0 , e y0 , e t0 , v x0 , u y0 . 

of the island resonances distort the tune dependence on the amplitude and 
the emittance growth caused by the island resonance. These effects are to 
some extent measureable, and can lead to a way to correct the resonance. 
Islands in 4-dimensional phase space, unlike islands in 2-dimensional phase 
space, are not easy to visualize. This study will show how to visualize the 
islands in 4-dimensional phase space by studying the tune dependence on the 
amplitude of the particle. Again, this effect is to some extent measureable, 
and can lead to a way to measure the strength of the resonance. 

Resonances in 2-dimenional phase space 

This section will study an island resonance in 2-dimensional phase space, the 
3v x = 19 resonance in the storage ring of the SNS (Spallation Neutron Source 
at Oak Ridge). This resonance is driven by the random error sextupole, 62 > 
in the magnets. In the absence of space charge forces, a non-linear field that 
produces an appreciable tune dependence on the amplitude of the particle 
motion is provided by the fringe field of the quadrupoles [T|, which causes 
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nuxvs. epxo, 3nux=1 9 resonance 
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Figure 2: v x versus e x o for the 3u x = 19 resonance excited by a ran- 
dom &2 in the SNS magnets. The initial coordinates have xq = 24mm 
to x = 44mm in steps of 2 mm with p x0 = and e y0 = .u x0 = 
6.3267,^0 = 6.2267. In the figure x0,px0, yO, pyO, epxO, eptO, nuxO, nuyO 
represent x , p x0 , y , p y0 , e x0 , e y0 , e t0 , v x0 , v y0 . 

depx.max vs. epxo, 3nux=19 resonance 
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Figure 3: de x>max vs e x0 for the 3v x = 19 resonance excited by a ran- 
dom b 2 in the SNS magnets. The initial coordinates x ,p x0 lie along the 
direction in Xq,p x o phase space given by p X Q = and e y = 0. u x0 = 
6.3267, u y0 = 6.2267. In the figure x0,px0, y0,py0, epxO, eptO, nuxO, nuyO 
represent x ,p x0 , yo,P y o, e^o, e y o, e y o, Qo, ^o, Vyo- 
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depx.max vs nux, 3nux=19 resonance 
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Figure 4: de x>max vs v x for the 3v x = 19 resonance excited by a random 
b 2 in the SNS magnets. The initial coordinates Xq,p x q he along the di- 
rection in Xo,Pxo phase space given by p x o = and e y0 = 0. u x q = 
6.3267, Vyo = 6.2267. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 repre- 
sent Xq, P X Q, HQ, PyQ, CrzO; tyQ, 6*0) V x0, ^y0- 

the tune to increase with the amplitude . Starting with a zero amplitude 
tune of v x = 6.3267 , 3u x = 19 — .02, one finds the phase space plot shown 
in Fig. ^ This plot was generated starting with xq = 24mm, p xo = 
and then increasing xq in steps of 2 mm keeping p xo = and tracking the 
particle for 1000 turns for each xq. The tune starts below the resonance, 
and increases with amplitude. The islands appear when the tune reaches 
v x = 6.3333. The islands are reached when the initial amplitude is increased 
to Xq = 30mm, p xo = 0, which corresponds to an initial horizontal emittance 
of about e X Q = 115 mm.mrad. 

The distortion of the tune dependence on amplitude is shown in Fig. |21 
where the horizontal tune is plotted against e^n- Fig.EJis generated by doing 
a search along the p xo = direction, increasing xq in steps of 2 mm. Fig. |21 
shows a flat region where v x = 6.3333 which will be seen to indicate the 
region where x ,p xo is crossing the island. When x ,p XQ are inside the island, 
the paticle motion is doing a slow oscilation around the fixed point at the 
center of the island. The particle motion then contains more than one tune, 
and the tune with the largest ampltude is the tune of the fixed point, 6.3333. 
Inside the island, the tune shown in Fig. El is this tune with the largest 
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amplitude. One can verify that the beginning and the end of the flat region 
occur at the same e X Q which correpond to the borders of the islands along the 
Px = direction. Fig. |21 can be used to find the width of the island along 
the p Xo = direction which is given by the width of the flat region in Fig. |21 
The islands generated by the 3z/ x = 19 resonance indicate a growth in 
the particle emittance. If one starts a particle inside the islands with the 
initial emittance e x o, then as the particle moves around the fixed point, its 
emittance will change and reach the maximun value of e XjTnax . One can use 
as a measure of the emittance growth the quantity 



• de x>max is plotted against e x0 in Fig. El It shows a maximun value for de x>max 
which can be used as a measure of the strength of the resonance. 

The results shown in Fig. El and Fig. |3] can to some extent be measured 
experimentally to find a way to correct this resonance with sextupole cor- 
rectors. An interesting plot that may corrrespond more closely to something 
that might be measured is to plot the emittance growth, de xmax vs the tune, 
v x . This plot can be found by combining Fig. [2] and Fig. Eland is shown in 
Fig. HI In this plot, there is a peak in the emittance growth which occurs 
when z/ x =6.3333, the resonance tune. The peak could be used as a measure 
of the resonance strength to correct the resonance. 

Using the flat region in the tune dependence on amplitude along the 
p x o = direction to indicate the strength of the resonance can sometimes 
lead to an error in setting the strengths of the sextupole correctors. Changing 
the strength of the correctors can cause the islands to move, and the width 
of the island as indicated by the flat region in the tune dependence on the 
amplitude may, for example, appear smaller because one is now crossing the 
island at a place where it is narrower. Also, when the islands move the the 
emittance growth , de Xymax may seem smaller. In order to avoid this error 
due to the movement of the islands when the corrector strengths are changed, 
one has to measure the tune dependence on amplitude along enough different 
directions in phase space so that one finds the largest width of the island and 
the largest emittance growth for the Xq,p x q inside the island. The same sort 
of argument also applies in finding the emittance growth dependence on the 
amplitude. In Fig. and Fig. |H1 results are given found by going out along 
many directions in phase space. One needs just enough directions to cover 
one of the three islands in Fig. ^ Here 6 directions were used in the region 
of phase space covered by one island. 
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Figure 5: v x vs e X Q for the 3u x = 19 resonance excited by a ran- 
dom b2 in the SNS magnets. The initial coordinates xo,p x o lie 
along 6 directions in the xo,p x o phase space and e y o = .v x0 = 
6.3267, Vyo = 6.2267. In the figure xO,pxO,yO,pyO,epxO,eptO,nuxO,nuyO repre- 
sent x , p x0 , yo, p y o, e x0 , e y0 , e t0 , v x(h v y o- 
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Figure 6: de x>max vs e x0 for the 3z/ x = 19 resonance excited by a 
random 6 2 in the SNS magnets. The initial coordinates Xo,p x o lie 
along 6 directions in the Xq,p x q phase space and e y o = 0. u x0 = 
6.3267, u y o = 6.2267. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 repre- 
sent Xo,p x Q, Uo,PyO, £x0, £y0, QO5 ^xO, v y0- 



nux vs. epxo, 3nux=19 resonance 
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depx.max vs nux, 3nux=19 resonance 

1.4 
1.2 
1 

X 

E 0.8 

X 

Q. 0.6 

CD 

~o 

0.4 
0.2 


6.328 6.33 6.332 6.334 6.336 6.338 6.34 6.342 6.344 6.346 
nux 

Figure 7: de x ^ max vs v x for the "iv x = 19 resonance excited by a 
random b 2 in the SNS magnets. The initial coordinates Xq,p x q lie 
along 6 directions in the xo,p x o phase space and e y o = 0. v x q = 
6.3267, UyQ = 6.2267. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 repre- 
sent xo,p x o, yoiPyo, e x o, e y o, eto, v x o-, v y o- 

In Fig. that plots v x versus e x0 , for each e x0 there are now 6 points plotted 
corresponding to the 6 directions. Points that have the tune, nux=6.3333, 
lie inside the island. The width of the island in e x o as given by the width of 
the points that lie on the nux=6.3333 line and is much larger than the width 
found using just the p x o=0 direction as it includes the direction that crosses 
the island close to where the island is widest. In the same way, Fig. |U]now 
shows the largest emittance growth for the xo,p x o that are inside the island. 

In applying the above results to develope a procedure for using the cor- 
rectors to correct a resonance, one will probably use a sample of the x ,p x0 
that is convenient for a particular ring and its injection system. One would 
then measure either the tune or the emittance growth of enough particles in 
the sample to be able to find the width of the island from the tune measure- 
ment or the largest emittance growth for paticles in the sample. One has to 
have a large enough sample that one is not misled by the movement of the 
islands when the excitation of the correctors are changed. To measure the 
largest emittance growth for particles in the sample, one can measure the fre- 
quency spectrum of the betatron oscillations. Particles inside an island will 
have the frequency that corresponds to the resonant tune 6.3333. The largest 
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Figure 8: u x vs e x o for the 3v x = 19 resonance excited by a random 62 in 
the SNS magnets and corrected using two sextupole correctors. The initial 
coordinates xq, p x o lie along 6 directions in the xq, p x o phase space and e y o = 0. 
v x q = 6.3267, Vyo = 6.2267. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 
represent x , p x o, Vo, P y o, e*o, <%o, %), v x o, Vyo- 



amplitude found in the frequency spectrum for the resonant frequency, can 
be used as a measure of the largest emittance growth. Something like this 
was done at RHIC 2J. The measurements can be done for a weak beam (no 
space charge effects) as it has been shown |3| that the resonance correction 
for resonances generated by magnet errors in the absence of space charge will 
also work fairly well in the presence of space charge. 

An interesting plot that may corrrespond more closely to something that 
might be measured is to plot the emittance growth, de Xjmax vs the tune, v x . 
This plot can be found by combining Fig. and Fig. Eland is shown in Fig. [7| 

The correction of the island resonances for the 3i/ x = 19 example can be 
done with two sextupole correctors properly located around the ring. The two 
correctors can then be adjusted one at a time, and for each setting of the 
correctors, the island width can measured from the dependnce of the tune 
on the amplitude as shown in Fig. El or the dependence of the amplitude 
growth on the tune as shown by combining Fig. and Fig. El Simulating 
this procedure gives the results shown in Fig. |H1 and Fig. El The emittance 
growth ,the maximun depx,max, is reduced by almost a factor of 10 by the 
best setting of the correctors that was found. 
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depx.max vs. epxo, 3nux=19 resonance 
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Figure 9: de x>max vs e x o for the 3u x = 19 resonance excited by a random fe 2 in 
the SNS magnets and corrected using two sextupole correctors. The initial 
coordinates xq,p x q lie along 6 directions in the xo,p x o phase space and e y o = 
■ Vxo = 6.3267, Vyo = 6.2267. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 
represent x , p x0 , yo, p y o, e*o, e y0 , e t0} v x0 , v y o. 

The theoretical result for the width of the island indicates that the width 
can be made zero when the correctors are set so that <ii/3n=0 where dv^ is 
the stopband [3] of the 3u x = 19 resonance due to the field errors in the 
magnets that are driving this resonance. According to the theoretical result 
the required setting of the correctors does not depend on the field producing 
the dependence of the tune on amplitude. This result is not useful for setting 
the correctors as the errors in the magnets are not known. However, this 
result helps to explain why the setting of the correctors in the absence of 
space charge fields will also correct the resonance when the space charge 
fields are present. It was found for the 3v x = 19 resonance that setting of 
the correctors to make di> 30 =0 will only reduce the emittance growth, the 
maximun de x>maxi by a factor of about 3 instead of the factor of 10 found 
above. It was found that for the SNS j2] this weaker correction that makes 
^30=0 was good enough to reduce the beam losses due to the resonance and 
space charge fields. 




9 



0.5 
0.45 

0.4 
0.35 
- 0.3 



0.2 
0.15 

0.1 
0.05 



Figure 10: The emittance spread de t vs e t o for the v x + 2u y = 19 
resonance excited by a random 62 in the SNS magnets . The initial 
coordinates x ,p x0 ,y ,Pyo lie along the direction in x ,p x0 , yo,P y o phase 
space given by p x0 = 0, p y0 = and e y0 = e x0 M x0 = 6.3933, 
UyQ = 6.2933. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 represent 

x 0,Px0, Vo,PyO, e x0, £y0, e i(b v x§i ^y0- 

Resonances in 4-dimenional phase space 

This section will study an island resonance in 4-dimensional phase space, the 
v x + 1v y = 19 in the storage ring of the SNS (Spallation Neutron Source at 
Oak Ridge). This resonance is driven by the random error sextupole, 62, in 
the magnets. Unlike, the 3u x = 19 resonance in 2-dimensional phase space, it 
is difficult to visualise the islands associated with the v x + 2v y = 19 resonance 
by looking at the particle motion in 4-dimensional phase space. Instead, we 
will look at the distortion due to the resonance on the dependence of the 
tune and the emittance growth on the amplitude of the particle motion. In 
Fig. [TUJ the spread in total emittance , 

)/(* 

.max ,'mvn) 

is plotted against the initial total emittance, e t0 . e tjmax is the largest total 
emittance and €t, m in is the smallest total emittance reached by particle with 
the initial coordinates, Xo,p x o, yo,Pyo, and initial total emittance, e t o,tracked 
for 1000 turns. Points were found along a particular direction in phase space 
given by p x o = 0, p y o = and e x = e y .The emittance spread becomes large 
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in the region e t0 =125 to e t0 =300, and this region may be used as a measure 
of the width of the island along this particular direction , p x0 = 0, p y0 = 
and e x = e y . The low value of de t near the center of region shows that this 
particular direction goes close to the center of the island where we expect 
de t =0. Fig, fm shows the emittance growth, 

det^max — (tt,max ~ Qo) / ' e t0 

as a function of e t Q. The goal of a correction scheme would be to reduce 
the maximun value of de t ^ m ax i n ths plot. In the 2-dimensional case of the 
3u x = 19 resonance, the plot of v x vs. e X Q had a flat region where u x had the 
constant value of 6.3333 which gave the width of the island. It was found 
for the u x + 2v y = 19 resonance, a flat region would be found if we plotted 
v x + 2v y vs. tto and the constant value in the flat region is v x + 2v y = 19. 
This is shown in Fig. The width of the flat region coincides with the 
width of the island as found from Fig. E3 which plots the emittance spread 
vs. €to- For coupled motion, the particle motion contains more than one 
tune, and what is plotted in Fig. is the main tune or the tune with the 
largest amplitude. 

An interesting plot that may corrrespond more closely to something that 
might be measured is to plot the emittance growth, de timax vs v x + 2v y . 
This plot can be found by combining Fig. ^2 and Fig. and is shown in 
Fig. In this plot, there is a peak in the emittance growth which occurs 
when v x + 2v y = 19. The peak could be used as a measure of the resonance 
strength to correct the resonance. 

As in the 2-dimesional case, to correct the v x + 2v y = 19 resonance , one 
has to repeat the above calculations for all directions in phase space that 
pass through a particular island, in order not to be misled by the movement 
of the islands when the correctors are changed. These results are shown in 
Fig- EJ Fig. IToland Fig. 1161 where for each e t o runs are done for 6 directions 
in Xq,p x o and for 6 directions in yo,p y o- In these figures there are, then, 36 
points for each value of e t o- The results indicated by these figures for the 
width of the island and the maximun emittance growth are not too different 
from those found above where only one direction in phase space was used, as 
the particular direction used appears to pass close to the center of the island. 

An interesting plot that may corrrespond more closely to something that 
might be measured is to plot the emittance growth , de tjma x vs u x + 2v y . 
This plot can be found by combining Fig. and Fig. and is shown in 
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Figure 11: The emittance sgrowth det^max vs eto for the v x + 2u y = 19 
resonance excited by a random b 2 in the SNS magnets . The initial 
coordinates xo,p x o,yo,p y o he along the direction in xo,p x o,yo,Pyo phase 
space given by p x0 = 0, p y0 = and e y0 = e x0 .u x0 = 6.3933, 
UyQ = 6.2933. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 represent 

x 0,Px0, U0,Py0, e x0, £y0, Qtb v x§i v y0- 
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Figure 12: v x + 2v y vs e t o for the v x + 2u y = 19 resonance ex- 
cited by a random b 2 in the SNS magnets . The initial coor- 
dinates xo,p x o,yo,Pyo lie along the direction in xo,p x o,yo,p y o phase 
space given by p x0 = 0, p y0 = and e y0 = e x0 .u x0 = 6.3933, 
u y Q = 6.2933. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 represent 

Xq,Px0, U0,Py0, Crf), e y0, Qo, ^xO, V y Q. 
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deptx.max vs nux, nux+2nuy=19 resonance 
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Figure 13: de t ^ max vs v x + 2v y for the v x + 2v y = 19 resonance 
excited by a random 62 in the SNS magnets . The initial co- 
ordinates x ,p x o,y ,p y0 lie along the direction in x ,p X Q,y ,p y0 phase 
space given by p x0 = 0, p y0 = and e y0 = e x0 .u x0 = 6.3933, 
VyQ = 6.2933. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 represent 

x 0,Px0, V0,Py0, e x0, £y0, e i(b v x§i ^y0- 

Fig. El In this plot, there is a peak in the emittance growth which occurs 
when v x + 2u y = 19. The peak could be used as a measure of the resonance 
strength to correct the resonance. 

The correction of the island resonances for the u x + 2i> y = 19 resonance 
can be done with two sextupole correctors properly located around the ring. 
Assuming that one has a way to measure either the width of the island or 
the maximun emittance growth using the above results as a guide , then the 
two correctors can be adjusted one at a time to reduce the emittance growth 
or the island width. Simulating this procedure gives the results shown in 
Fig. [__ and Fig. [THl The emittance growth ,the maximun de tt max, is reduced 
by almost a factor of 15 by the best setting of the correctors that was found. 
In this case, setting the correctors to zero the stopband, dui2=0, reduces the 
maximun emittance growth by a factor of 5. 

Higher order resonances 

In the above, results were presented for the 3nux=19 and the v x + 2u y = 19 
resonances. Similar studies have also been done for all four of the third 
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Figure 14: The emittance spread de t vs e t0 for the v x + 2v y = 
19 resonance excited by a random 62 in the SNS magnets . The 
initial coordinates xo,p x o,yo,Pyo he along 6 directions in the xo,p x o 
phase space and 6 directions in yo,p y o space . u M = 6.3933, 
UyQ = 6.2933. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 represent 

x 0i PxO, V0, PyO, e x0, £y0, e t0, v x§i ^y0- 
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dept,max vs. epto, nux+2nuy=19 resonance 
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Figure 15: The emittance sgrowth det, m ax vs e t o for the v x + 2u y = 
19 resonance excited by a random hi in the SNS magnets . he 
initial coordinates xo,p x o,yo,P y o he along 6 directions in the xo,p x o 
phase space and 6 directions in yo,p y o space . u x0 = 6.3933, 
UyQ = 6.2933. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 represent 

x 0iPx0i VdiPyQi e x0, £y0, e t0i v x§i ^y0- 
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Figure 16: v x + 2v y vs e t0 for the v x + 2v y = 19 resonance excited by a 
random 6 2 in the SNS magnets . The initial coordinates Xo,p x o,yo,p y o lie 
along 6 directions in the xo,p x o phase space and 6 directions in y0,py0 space 
■ v x o = 6.3933,1/^0 = 6.2933. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 
represent x Q , p x0 , yo, p y o, e^o, e y o, Qo, "xo, Vyo- 
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deptx.max vs nux, nux+2nuy=19 resonance 
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Figure 17: de t ,max vs ^ + 2z/ y for the v x + 2^ = 19 resonance excited by 
a random 62 in the SNS magnets . The initial coordinates XQ,p x o,yo,P y o lie 
along 6 directions in the xo,p x o phase space and 6 directions in yo,p y o space 
■ Vxo = 6.3933,1^0 = 6.2933. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 
represent x , p x0 , y , p y0 , e x0 , e y0 , e t0 , v x0 , u y0 . 
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Figure 18: The emittance sgrowth de ty max vs e t o for the v x + 2v y = 19 res- 
onance excited by a random 62 in the SNS magnets and corrected using 
two sextupole correctors. The initial coordinates xo,p x o,yo,P y o he along 
6 directions in the xO, pxO phase space and 6 directions in yo,P y o space 
■v x o = 6.3933,1/^0 = 6.2933. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 
represent x , p x0 , y , P y o, e x0 , e y0 , e t0 , u x0 , u y0 . 
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order resonances mv x + nv y = 19, m + n = 3, and for all five of the fourth 
order resonances mv x + nv y = 25, m + n = 4, with similar results. It is our 
expectation that similar results would be found for the higher order resonance 
mu x + nu y = q, where m,n,q are integers. In particular, if mu x + nu y is plotted 
versus e t o, this plot will contain a flat region where mu x + nu y is costant at 
mv x + nu y = q, and the width of this flat region will give the width of the 
islands associated with the mu x + nu y = q resonance. 
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nux+2nuy vs. epto, nux+2nuy=19 resonance 



=> 19.01 
c 

+ 19.005 

X 

=> 19 
c 

18.995 
18.99 
18.985 















i - 












! 










( 












1 

1 














• 
























- 



























50 100 150 200 250 300 350 400 450 
eptO mm.mrad. 



Figure 19: v x + 2\nu y vs e t0 for the v x + 2v y = 19 resonance ex- 
cited by a random b 2 in the SNS magnets and corrected using two sex- 
tupole correctors. The initial coordinates x ,p x0 ,y ,p y0 lie along 6 direc- 
tions in the xo,p x o phase space and 6 directions in yo,p y o space .v x o = 
6.3933,^0 = 6.2933. In the figure X0,px0,y0,py0,epx0,ept0,nux0,nuy0 rep- 
resent Xo,p x o, l/0,Py0, CxO, £y0, CtO, V x0 , U y0 . 
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